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We report on a technique to improve the continuous loading of atomic strontium into a magnetic trap from a 
Magneto-Optical Trap (MOT). This is achieved by adding a depumping laser tuned to the —» ^5 1 (688-nm) 

transition. The depumping laser increases atom number in the magnetic trap and subsequent cooling stages by 
up to 65 % for the bosonic isotopes and up to 30 % for the fermionic isotope of strontium. We optimize this trap 
loading strategy with respect to the 688-nm laser detuning, intensity, and beam size. To understand the results, 
we develop a one-dimensional rate equation model of the system, which is in good agreement with the data. We 
discuss the use of other transitions in strontium for accelerated trap loading and the application of the technique 
to other alkaline-earth-like atoms. 


I. INTRODUCTION 

Alkaline-earth-like (AE) atoms have received a great deal 
of recent interest due to the distinctive properties of their level 
structure [1, 2]. The largely disconnected singlet and triplet 
states in these atoms give rise to very narrow optical transi¬ 
tions, which could form the basis for an improved time stan¬ 
dard [3, 4]. These transitions are also advantageous in a wide 
variety of other applications. For example, their low photon 
scattering rates allow for the production of highly-excited Ry¬ 
dberg atoms with reduced decoherence compared to alkali- 
metals [5]. Magnetic held insensitive singlet and triplet lev¬ 
els make AE atoms attractive for precision measurement and 
quantum sensing applications [6, 7]. In fermionic isotopes, 
these states manifest SU(2/ + 1) spin symmetry, where I is 
the nuclear angular momentum, allowing quantum simulation 
of hamiltonians that are inaccessible with alkali atoms [8-10]. 
All of these applications require or beneht from a combination 
of large atom number and short experimental cycle times. 

Recent advances in cooling and trapping techniques en¬ 
abled production of the hrst strontium degenerate gases [11- 
15]. The small negative s-wave scattering length of the most 
abundant isotope, ^®Sr, hampered initial efforts to create Bose- 
Einstein condensates [16-20]. While the other stable isotopes 
(87sr, 86 §j.^ ®"^Sr) possess favorable scattering lengths, 

their low natural abundance initially prevented Magneto- 
Optical Traps (MOT) from collecting enough atoms to reach 
degeneracy. Fortuitously, laser cooling of strontium on the 
461-nm line populates a magnetically-conhned, metastable 
reservoir of atoms in the ^P 2 state (see Figure 1) [21]. The 
long lifetime of this reservoir (typically >10 s) compared to 
the MOT allows for the accumulation of sufficient populations 
of ®’Sr, ®®Sr, or ^'^Sr for forced evaporation or sympathetic 
cooling of *^Sr [11-15]. The pK temperatures attainable 
with laser cooling on the 689-nm, intercombination transition 
(Figure 1) lead to short evaporation times to reach degener¬ 
acy. Given the low abundance of the interacting isotopes, this 
means that the reservoir loading time usually dominates the 
experimental cycle [11, 13, 14, 22]. 

Typical Sr degenerate gas experiments hrst use a MOT op¬ 
erating on the *5o —> *Ri, 461-nm transition (bMOT) to cap- 
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FIG. 1. (Color online) The low-lying energy levels of bosonic stron¬ 
tium with linewidths taken from [1, 23, 24]. A Magneto-Optical Trap 
operating on the blue, 461-nm transition (bMOT) captures atoms 
from a Zeeman-slowed beam. Atoms in the hMOT continuously leak 
into the long-lived ^P 2 state, which is magnetically trapped by the 
hMOT quadrupole field. Two lasers at 688 nm and 679 nm increase 
the magnetic trap loading rate by pumping atoms that populate ^Pi 
into ^P 2 . The 679-nm laser and a 707-nm laser return atoms to the 
ground state via the state once magnetic trap loading is com¬ 
plete. A Magneto-Optical Trap operating on the red, 689-nm transi¬ 
tion (rMOT) then cools the sample to s;! pK. 


ture atoms from a Zeeman-slowed atomic beam and cool them 
to 1 mK. Atoms slowly leak out of the bMOT cooling cycle 
(1:50,000 branching ratio) and into the metastable mani¬ 
fold, where they populate the ^P 2 and ^Pi states in a 1:2 ra¬ 
tio [23, 24]. The bMOT quadrupole held can magnetically 
trap atoms in the state (the Lande g-factor gj = 3/2 for 
bosonic isotopes, where J is the electronic angular momen¬ 
tum). Repumping lasers return ^P 2 atoms to the ground state 
once magnetic trap loading is complete, which, depending on 
the isotope, can take 30 s or more [11, 13, 22]. Loading times 
are also long for experiments with isotopic mixtures, since 
the isotope shifts of the 461-nm transition are on the same 
order of magnitude as the linewidth [14, 22]. This prohibits 
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FIG. 2. (Color online) Atom number enhancement as a function of 688-nm laser detuning, Ajgg, for all strontium isotopes. The detuning zero is 
referenced to **Sr. For the data shown, the rMOT recaptures the less abundant isotopes and the hMOT recaptures *®Sr. The depumper saturation 
parameter was isgg « 35 (~ 50) for **Sr Sr,*'* Sr) and the 679-nm laser detuning, A(,jg, was set to maximize hMOT fluorescence. 

We label the fermionic hyperfine transitions with F^F’, where F is the total angular momentum quantum number for ^Pi and F' the 
corresponding quantum number for ^5 1 . Inset: The detuning-dependent scattering rate for each transition between and i Zeeman levels, 
averaged over the volume of a one-dimensional hMOT with iggg = 1 (see Section IV). Solid blue curves pump to a ^5 j Zeeman level that 
can decay to a magnetically trappable ^f ’2 Zeeman state, but dashed red curves do not. Asymmetric lineshapes arise in the atom number 
enhancement because the dashed red scattering rate curves dominate at negative detuning. 


efficient simultaneous loading of the magnetic trap. A sec¬ 
ond stage Magneto-Optical Trap using the 689-nm, intercom¬ 
bination line (rMOT) cools these atoms to =» 1 fjK and facili¬ 
tates loading into an optical dipole trap. Evaporation proceeds 
quickly due to the low initial temperature and degeneracy can 
be reached in :«1 s for most isotopes [22]. 

Here we present a technique to reduce the reservoir loading 
time or, equivalently, increase the atom number for experi¬ 
ments with strontium as first suggested in [25]. The method 
relies on continuous optical pumping of atoms from the short¬ 
lived state into the magnetically trapped reservoir us¬ 
ing the —» ^5 1 , 688-nm transition. This greatly reduces 

the steady-state atom number in the hMOT, but increases the 
flux of low-field seeking atoms into the metastable reservoir. 
Although the ^P 2 '.^P\ branching ratio from the *^2 state sug¬ 
gests that atom number should be enhanced by a factor of 
three (see Figure 1), we show that this estimate is incorrect 
since it does not consider the reduction in hMOT atom num¬ 
ber caused by the 688-nm laser. 

We describe our experimental apparatus in Section 11 with 
an emphasis on the details relevant for the accelerated load¬ 
ing scheme. Section 111 explains the measurement procedure 
and results. In Section IV, we develop a rate equation model 
and demonstrate that our data is in agreement with expecta¬ 
tions. We also simulate the trap loading enhancement for sev¬ 
eral other transitions in strontium and two in calcium. Section 
V is a summary of our results and outlook. 


II. APPARATUS 

Our experimental setup is similar to other strontium appa¬ 
ratuses designed for optical clock and degenerate gas exper¬ 
iments [25-28]. An oven with a microtubule array nozzle, 
heated to 600 °C, creates an atomic strontium beam. Two 
stages of differential pumping prevent the pressure in the ex¬ 
periment chamber (6 x 10“^' Torr) from rising while the oven 
is in operation. The atomic beam passes through a transverse 
cooling stage, which consists of two orthogonal, retroreflected 
461-nm laser beams. Each beam has sslO mW of power, a 1:3 
aspect ratio (1 /e^ radius of 9 mm along the atomic beam axis), 
and -10 MHz detuning from the *5o ^ ^Pi transition. The 
Zeeman slower is a 35-cm long, multilayer, variable pitch coil 
located immediately after the transverse cooling stage. It is 
pumped with a!48 mW of -600 MHz detuned 461 -nm light, 
which is focused onto the oven nozzle with an initial 1/e^ ra¬ 
dius of 5 mm. 

The hMOT has a standard retroreflected, six beam config¬ 
uration. Each beam has a 1/e^ radius of 8 mm, a detuning 
A 461 = -45 MHz, and contains either ^7 mW (for bosonic 
isotope data) or ^9 mW (for ^^Sr data) of power. These pa¬ 
rameters give i 46 i = I/Isat ~ 0.16 per beam for the bosons 
and i 46 i ~ 0.21 per beam for the fermion. The quadrupole 
coil has a vertical axis of symmetry and produces a magnetic 
field gradient of 6 mT/cm along that axis during hMOT oper¬ 
ation. The hMOT field gradient is sufficient to magnetically 
trap ^P 2 atoms in the low-field seeking \mj = 1) and \mj -2) 
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FIG. 3. (Color online) N/Nq for **Sr recaptured in the bMOT as a 
function of A 579 . The 688 -nm laser has Jsss ~ 35 and Aggg 30 MHz. 
We reference the 679-nm detuning to the bMOT fluorescence max¬ 
imum, indicated by the dashed vertical line. As shown in the inset 
to Figure 2, the asymmetry arises from detuning-dependent scatter¬ 
ing rates. Transitions to 1 Zeeman levels that can decay into the 
magnetic trap are predominately blue detuned, whereas red detuned 
transitions populate levels that decay to high-held seeking ^P 2 Zee- 
man states. 


Zeeman sublevels. In our vacuum chamber, the position of 
two recessed viewports along the symmetry axis of the coils 
limits the trap depth for the \mj - 1) state to xi5 mK. This 
limitation is unimportant for us since our bMOT loads ssl mK 
atoms into the magnetic trap, but it suggests that experiments 
with larger vacuum chambers may find that a higher tempera¬ 
ture bMOT optimizes magnetic trap loading [25]. 

Two repumping lasers addressing the 679-nm, ^Pq 1 
and the 707-nm, ^P 2 1 transitions are used to return 

atoms in the magnetic trap to the ground state. The two beams 
co-propagate with the Zeeman slower beam, share a 1/e^ ra¬ 
dius of =sl cm, and contain =s2.5 mW (679 nm) and =s4.5 mW 
(707 nm) of power. For experiments with the bosonic iso¬ 
topes, we lock the repump laser frequencies using slow feed¬ 
back from a HighFinesse WS7 wavemeter [29]. The locking 
stability is +5 MHz, which is much narrower than the ob¬ 
served bosonic repumping linewidth [30]. The presence of 
hyperfine structure in the fermion complicates repumping on 
the ^P 2 —> 1 transition. In order to cover as much of the 

=s5.5 GHz hyperfine spectrum of the transition as possible, we 
modulate the 707-nm laser frequency at ^700 Hz. To further 
increase coverage of the hyperfine spectrum, we use a second 
707-nm laser that we modulate at ss600 Hz. When optimized, 
application of the second laser to the experiment increases the 
®^Sr atom number by about 10 %. For the fermionic data, the 
679-nm laser is locked to the \^Po,F - 9/2) —> = 

11 /2) transition (where F - I + J) using the wavemeter. 

The linewidth of our 689-nm master oscillator is stabilized 
below the natural linewidth of the —» ^P\ resonance us¬ 


ing a Pound-Drever-Hall lock to an optical cavity (finesse 
ai240,000) [31]. We injection lock a slave laser diode to the 
master to obtain sufficient power for trapping (for the fermion, 
this laser pumps the \ ^So,F - 9/2) ^ = 11/2) transi¬ 

tion). Dichroic beamsplitters overlap the 689-nm light for the 
rMOT with the bMOT beams. In each rMOT arm, the power 
is ss3.5 mW and the 1/e^ radius is 2.5 mm. The quadrupole 
field gradient switches to 0.16 mT/cm for rMOT operation. 
For the first 100 ms of rMOT operation, we frequency mod¬ 
ulate the trapping laser at 30 kHz with a modulation depth 
of 1 MHz to increase the capture velocity of the rMOT. Over 
the next 400 ms, we linearly reduce the modulation depth to 
100 kHz while simultaneously ramping the optical power to 
100 fiW with a half-gaussian temporal profile. In this work, 
we terminate rMOT operation at this stage (T « 2 fiK), but we 
can cool further by turning off the frequency modulation and 
reducing the intensity. The hyperfine structure of the fermion 
means that a second slave laser is required to make a stable 
rMOT [32]. A beatnote lock to the master stabilizes the sec¬ 
ond laser to the = 9/2) ^ pFi,F' = 9/2) line [33]. 

The second slave laser provides ^ 800 /iW of light per rMOT 
beam to the experiment. Aside from the reduction in ini¬ 
tial power, the intensity and modulation ramps are identical 
to those for the trapping laser. We also linearly increase the 
rMOT magnetic field gradient to 0.24 mT/cm during the final 
400 ms of the rMOT when trapping *^Sr to increase the atomic 
density. 

To enhance magnetic trap loading, a 688-nm laser resonant 
with the ^Pi —> 1 transition pumps atoms that decay to ^Pi 
out of the bMOT cycle and into We call this laser the 
depumper because it makes the bMOT transition less closed. 
The depumper is a Littman-Metcalf configuration laser which 
we built using a laser diode (model HL6738MG [29]) that was 
AR-coated in-house [34]. The laser provides up to 2.8 mW of 
light to the experiment, which corresponds to s(,ss - 11 hat ~ 
50. The 688-nm beam enters the chamber horizontally and 
perpendicular to the Zeeman slower axis. It has a 1/e^ radius 
W688 = 1 -35 mm except where otherwise noted. We stabilize 
the 688-nm laser detuning, A688, to within +3 MHz by locking 
to the wavemeter. 

We measure the magnetic trap loading enhancement by 
interleaving shots with the 688-nm laser on and off. The 
679-nm repumping laser closes the \ —» ^Pq leak to in¬ 
crease the depumper’s effect. It remains on during both shots 
of a depumper on/off pair of experimental runs, but does not 
affect trap loading when the 688-nm laser is off since the 
bMOT does not populate ^Pq. After 0.5 s to 30 s of reservoir 
loading, an acousto-optic modulator extinguishes the 688-nm 
beam and optical shutters open to allow 707-nm light to reach 
the experiment. Either the bMOT or the rMOT can recap¬ 
ture atoms from the magnetic trap for detection and imaging. 
However, rMOT recapture greatly improves signal-to-noise 
for ®^Sr, ^®Sr, and ^'^Sr, so we use it exclusively for these iso¬ 
topes. The bMOT recapture stage lasts 100 ms and rMOT re¬ 
capture consists of the full rMOT cycle described above. We 
take an absorption image after a 1 ms (25 ms) ballistic ex¬ 
pansion for bMOT (rMOT) recapture using a resonant, 10 //s 
pulse of 461-nm light with I j hat ~ 0.04. Numerical Integra- 
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FIG. 4. (Color online) N/Nq for **Sr recaptured in the bMOT ver¬ 
sus jggg and W 588 , where Wggs is the 1 /e^ radius of the depump beam. 
The standard errors are omitted for clarity, but are < 0.03. The op¬ 
timal trap enhancement occurs when >^688 roughly matches the l/e 
radius of the bMOT and ^688 ~ 1. ^688 and A(,^g are both sjSO MHz 
( corresponding to the maxima in Figures 2 and 3). 


tion of the image yields the atom number for each shot. 


III. RESULTS 

We study the magnetic trap loading enhancement as a func¬ 
tion of isotope, power, detuning, and beam size. The enhance¬ 
ment is measured by comparing the atom number recaptured 
in the rMOT or bMOT with and without the 688 -nm laser. 
We hnd that the depumper’s effect is independent of which 
MOT we use for atom recapture. The magnetic trap loading 
enhancement is given by the normalized atom number, N/Nq, 
where N is the atom number with the depumper on and No the 
number with it off. 

We investigate the loading enhancement as we scan the 
depumper across the 688 -nm transition. For this data set, 
we set 5688 ~ 35 for ^^Sr and 5688 ~ 50 for all other iso¬ 
topes. The repump laser frequencies are locked to maximize 
bMOT fluorescence. The magnetic trap loading time, fioad, 
for {*^^Sr, *^Sr, ^®Sr, ^^Sr) is {1.5 s, 10 s, 6 s, 7.5 s) result¬ 
ing in typical No of (2 x 10^,5 x 10®, 1 x 10’, 8 x 10®) in 
the rMOT. Adjustment of the wavemeter lockpoint allows us 
to scan A 688 - Figure 2 shows the depumping spectrum for 
all isotopes and hyperfine transitions, the locations of which 
are in good agreement with [35]. On average, we observe 
trap loading improvements of a!50 % for bosonic isotopes and 
ai25 % for the fermionic isotope even without detailed opti¬ 
mization of the depumping parameters. The peak enhance¬ 
ment for each isotope occurs when A 688 ~ 30 MHz from res¬ 
onance. In Figure 2 we also see that the choice of hyperfine 
transition is crucial for atom number gains in ^’Sr. Pump¬ 
ing to f'S i,F-7/2) and 1 ^ 51 ,^ = 9/2) is always detrimen¬ 
tal because these manifolds decay with >60 % probability to 
|®P 2 , F — 112) and \^P 2 , F = 9/2), which have Lande g-factors 
too small for magnetic trapping at the bMOT field gradient. 
Pumping to 11/2) yields a lineshape similar to that 

of bosonic isotopes, but with reduced amplitude. 

The asymmetric lineshapes observed in Figure 2 are due 
to the non-uniform magnetic fields in the bMOT. Quadrupole 


fields shift low-field seeking states to higher energy and high- 
field seeking states to lower energy. Because the Lande g- 
factor for 1 is larger than for ^Pi, this effect causes a 
blueshift for most transitions to 1 states that can decay to 
'^P 2 ,mj — 1) or '^P 2 ,mj - 2) (see inset to Figure 2). For the 
same reason, transitions to 1 states that can only populate 
untrapped Zeeman states are redshifted. As a result, trap 
loading is enhanced to the blue of resonance and reduced to 
the red of resonance. 

For this depumping scheme, application of the 679-nm laser 
during the bMOT is crucial; removing the 679-nm laser re¬ 
sults in ss50 % reduction of the enhancement. The effect of 
the 679-nm laser was studied by varying its detuning, A 679 . 
As shown in Figure 3, setting A 679 ~ 30 MHz (relative to 
the detuning that maximizes the bMOT fluorescence) adds an 
additional =sl0 % to the enhancement. The asymmetric line- 
shape is caused by the same mechanism discussed above for 
the 688 -nm transition. 

With A 688 and A 679 stabilized at their optimized values, we 
study the trap loading enhancement for *®Sr as a function 
of 5688 and W 688 - Slight focusing/defocusing of the 688 -nm 
beam changes the waist at the location of the bMOT, but the 
Rayleigh range is always larger than the bMOT l/e radius, 
^bMOT, for the parameter range we study. Figure 4 shows 
that trap loading enhancement increases with 5688 provided 
W 688 ^ ^'bMOT, with the optimal enhancement occurring when 
W 688 - ^bMOT- High 5688 increasingly reduces NjNo for larger 
beam waists. The data suggest that, for our bMOT parame¬ 
ters, a substantial number of atoms populate the manifold 
before being fully captured by the bMOT. These atoms exist 
outside the bMOT radius and are too hot for magnetic con¬ 
finement, but they are cold enough that they do not leave the 



FIG. 5. (Color online) Atom number enhancement versus 5588 with 
W688 = 1.35 mm. *’Sr data are for the 6"= 11/2 —» F' = 11/2 transi¬ 
tion and all isotopes are recaptured in the rMOT. Error bars represent 
the standard error in the mean for > 10 measurements. The 688-nm 
and 679-nm laser detunings are set to x: 30 MHz. 
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bMOT capture volume during the ms decay time for the 
—»'D 2 —»^Pi —»'5o path. The effect of varying sess and 
W 688 in the other isotopes was similar to the **Sr results. In 
Figure 5, we plot N/Nq for a wider range of the saturation pa¬ 
rameter at the optimum W 688 - All isotopes exhibit a steep rise 
in trap loading enhancement for iess ^ 1 , followed by a shal¬ 
low rolloff for S 688 > 1 ■ We find that the enhancement is very 
sensitive to the 688 -nm beam alignment and that the peak at 
■*688 ~ 1 is present only when the beam traverses the center of 
the bMOT. 

Before recapturing atoms from the magnetic trap, we do 
not first discard ground state atoms remaining in the bMOT. 
This increases both N and Nq, but decreases their ratio. This 
choice biases our results toward lower enhancement values, 
particularly for short load times and for ^*Sr. However, the re¬ 
duced N/No is the appropriate metric for evaluating the load¬ 
ing enhancement in most experiments, since the cycle time is 
typically limited by N. In situations where bMOT atoms are 
lost before recapturing from the magnetic trap, which arise in 
experiments with isotopic mixtures, the depumping technique 
is even more useful. If we remove the bMOT atoms before 
imaging, N/No increases by up to 15 %. 

Our technique reduces the trap loading time necessary to 
achieve a given atom number. For short trap loading times, 
N/No is a measure of the increased loading rate achieved with 
the depumping laser. This is the regime shown in Figure 2. 
For longer trap loading times, the atom number will saturate. 
Experiments requiring very large atom numbers close to the 
saturation limit can expect even greater reductions in loading 
time than suggested by the initial loading rate. We demon¬ 
strate this by fitting N(t) and No{t) with N(t) - 
where a is the loading time constant and is the sat¬ 
urated atom number (see inset to Figure 6 ). Inverting the 
fitted function yields the loading time necessary to reach a 
given atom number with the depumper on, t{N), or with the 
depumper off, fo(Ao). We plot the loading time reduction fac¬ 
tor, LTRF(A?/A?”“) = to(No^N)/t(N), for **Sr and ^^Sr in 
Figure 6 . The loading time reduction diverges as ^ 
since For example, to reach an atom number of 

the depumping technique can reduce the loading time 
by a factor of ss3. 


IV. SIMULATION 


To better understand the enhancement, we develop a one¬ 
dimensional rate equation model to simulate the bMOT 
depumping process. This was motivated by two features of 
our data: the asymmetric lineshapes depicted in Figure 2, and 
the discrepancy between the observed performance and the 3x 
initial estimate given by the ^P 2 '.^P\ branching ratio. A sim¬ 
ple calculation, based on analysis of the cascade of Clebsch- 
Gordan (CG) coefficients connecting 'Fi to ^Pi and ^P 2 , sug¬ 
gests that 3x trap loading enhancement is unlikely. How¬ 
ever, this CG calculation depends sensitively on the relative 
populations of the 'Fi Zeeman sublevels, which are position 
dependent, and the steady state atom number in the bMOT. 
Both of these complications prevent analysis of experimental 
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FIG. 6. (Color online) The Loading Time Reduction Factor (LTRF) 
for **Sr (dotted blue line) and ^''Sr (solid cyan line). The saturated 
atom number, is the asymptote of the trap loading curve with 
the 688-nm laser off (for **Sr, red triangles, inset). The loading time 
necessary to transfer N atoms into the rMOT with the depump laser 
on (off) is given by t (to). Inset: The raw data and fits associated 
with the ®®Sr (dotted blue) LTRF curve. The standard error is smaller 
than the data points. We find t and to by inverting the appropriate fit 
function (see text). 


performance by this method. A full simulation of the opti¬ 
cal pumping dynamics resolves both of these issues, allowing 
direct comparison of data with theory. Straightforward mod¬ 
ifications of the rate equation model allow us to compare our 
technique to alternative depumping transitions. 

In the rate equation model, we track the population, 
in each magnetic sublevel of i 6 {*5o, *Fi,... ,^5i), with m, 
the spin projection along the axis of a one-dimensional bMOT. 
Each level decays at a rate given by the appropriate linewidth, 
Yij, from Eigure 1 , 


Y<decay 

niy) 



( 1 ) 


where {i,tni-, l,my\j,mj) is the CG coefficient. In addition to 
the transitions shown in Eigure 1, we also include the ^D 2 —> 
'5o quadrupole decay because its linewidth is non-negligible 
compared to decay rates into the ^P states [23]. Since we are 
not interested in individual atom trajectories, we average the 
driven excitation rate for I/, m,) —> I /, m,), Ff.’^'^ , ,. ,, over 

the position and velocity distribution of the MOT, 
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FIG. 7. (Color online) The *®Sr curve from Figure 2 (blue cir¬ 
cles) plotted with simulation results (magenta line) from the one¬ 
dimensional rate equation model described in Section IV. Parame¬ 
ters for the simulation are identical to experimental conditions (see 
Section II). The model reproduces the qualitative lineshape and the 
peak enhancement of the data. The difference in the amplitude and 
width between the curves implies that our one-dimensional model 
does not accurately capture the full three-dimensional nature of our 
experiment. 

where Sjj is the saturation parameter, cr{my) is the fraction of 
Sij with polarization rriy e {-1,0,1), and the effective detun¬ 
ing between |i,m,) and |y,m^), includes Doppler and 

Zeeman shifts. We choose Xmax, Vmax to be much larger than 
the characteristic scale of p{x, v). The symmetry of a one¬ 
dimensional MOT permits us to model only the ;c > 0 region 
with all scattering rates then multiplied by two. This choice 
simplifies the tracking of magnetically trapped atoms because 
the magnetic field does not change sign in the simulation vol¬ 
ume. Taking the transformation —> 2s,y in Equation (3) 

while maintaining Em^cr(my) = 1 gives the correct scattering 
rate for the two bMOT beams. We assume a pure circular po¬ 
larization for both bMOT beams and a random polarization 
for the repumper and depumper. 

We describe the evolution of the populations, (Eiy,,, ), with a 
system of coupled differential equations, 

Pi,mi - J^'So^i,'So 

(rexc p _ r«c p 

Z_J \ ^ my)~^\j,mjy ‘Mi 

jMjMy 

r-idecay p j^decay p \ 

my) \i,mi)^\j,mj\ tUy) J ’ 

where is the Kronecker delta, 'Ti is the atomic flux from 
the Zeeman slower, and yS « 0.1 s ' is the experimentally mea¬ 
sured 1-body loss rate (the effect of which is negligible for 
states with short lifetimes). Without repumping, the bMOT 
loading time (< 100 ms) is short compared to the magnetic 
trap loading time, so we take = 0 for all i + ^P 2 . 

We solve algebraically for {P^Py- 2 , ■ ■ ■, P^P 2 , 2 ] numeri¬ 
cally integrate the resulting first-order equations from f = 0 
to t - fioad- The sum {Pip^a + P‘>P2,i + ^’So.o) gives the total 
atom number at f = fioad (the population of other states is neg¬ 


ligible), which we equate with N or No depending on whether 
the 688 -nm laser is on or off. A fit of the model to the ^*Sr, No 
versus hoad data, with s iggg = 0 and 'Fis„ as the only 

free parameter, matches the experiment to better than 4 % for 
all reservoir loading times (all other parameters are taken from 
Section II). We use the extracted value of TiSo for subse¬ 
quent simulations, but we find that the results are independent 
of !Fi So and p. 

We plot the simulated and measured NjNo for ®*Sr in Fig¬ 
ure 7. Parameters for the simulation are taken from Section 
II except for sggg and Agvg, which are the same as given in 
Figure 2. The simulation agrees reasonably well with experi¬ 
ment given the simplicity of the model and the absence of free 
parameters. It approximately reproduces the asymmetric line- 
shape and the magnitude of the peak trap loading enhance¬ 
ment. The one-dimensional model also qualitatively repli¬ 
cates the behavior of NjNo as a function of soss and Agvg. The 
difference in dimensionality between the ID simulation and 
3D experiment causes the mismatch in both the width and am¬ 
plitude of the lineshapes in Figure 7. The three-dimensional 
MOT beam configuration and magnetic quadrupole field com¬ 
plicate the optical pumping dynamics. 

The choice ofa7=l —» 7' = 1 transition as our depump¬ 
ing line potentially limits the trap loading enhancement, since 
the 688 -nm line has position-dependent dark states and small 
CG overlap with \^P 2 ,mj — 2). Furthermore, this transition 
requires a secondary laser to depopulate the ^Po state. Many 
repumping strategies exist for strontium and each of these pos¬ 
sesses a nearby depumping resonance [25, 36-38]. We as¬ 
sess the relative merit of the various schemes by simulating 
them with optimum parameters (Figure 8 ). The 5s5p ^Pi —> 
5s5d^D2 line at 487 nm and the 5s5p ^Pi —> 5s6d^D2 line 
at 397 nm have similar performance to the 688 -nm line. All 
other transitions for which linewidth data are available give 
less enhancement. For the SsSp^Pi —> 5p^^P2 transition, 
unfavorable relative Fande g-factors between the excited state 
and ^Pi marginally reduce the trap loading improvement. The 
linewidth of the —» 5 s 4 (i^D 2 transition is too narrow 

for efficient optical pumping at bMOT temperatures. 

We investigate the utility of the depumping scheme for 
other AE atoms. For Cd, Hg, Yb, Be, and Mg, the *D 2 state 
lies above the *Pi state, so efficient continuous loading of the 
metastable reservoir does not occur [2, 25]. Direct pumping to 
the magnetically trapped state is possible for these atoms [39]. 
The *D 2 : *5 0 branching ratio in Ba and Ra is very large, which 
means that cooling to temperatures below the magnetic trap 
depth may not be possible without repumping [40, 41]. The 
level structure of calcium combines several features that make 
depumping more effective than in strontium (see Figure 8 ). 
The ^P 2 '?Pi branching ratio is « 1:3 and, more importantly, 
the *D 2 ^ '5o quadrupole transition linewidth is compa¬ 
rable to 'D 2 —> ^Pj decay rates. Ca can be trapped in a 
MOT operating on ^P 2 ^Dj, transition [42], the loading 
of which could also benefit from this depumping technique. 
The loading enhancement for MOTs does not benefit from the 
detuning-dependent asymmetry seen for magnetic trap load¬ 
ing, which limits the simulated improvement in Ca to a!50 %. 
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FIG. 8 . (Color online) Simulations of several depumping schemes 
for bosonic strontium (solid lines) and calcium (dashed lines), which 
we label with the depumping transition excited state (the lower state 
is always 5s5p ^Pi for Sr and 4j4p for Ca). In all simulations, 
we use our ®*Sr values for rtMor, head, A 461 , and p. The 
depumper saturation parameter is = 1 except for the simula¬ 

tion of 5s4d^D2, where Sdepump = 2000 , which requires much higher 
saturation due to the narrow transition linewidth. For the simulation 
of SjOi" ^5 1 , A 579 is set to its optimal value. The Sr 5s5d ^£> 2 , Sr 
5s()d ^£> 2 , Sr 5p^ ^£* 2 , Ca Ap^ ^£ 2 , and Ca AsAd ^£>2 states may indi¬ 
rectly decay to ^Pq via intermediate states outside of the ^P manifold. 
The model ignores these processes, but atoms decaying into ^Pq can 
be recovered using e.g. a 679-nm laser (for Sr). The apparent offset 
of the enhancement for Ca is a gaussian pedestal with a full width at 
half maximum of approximately 500 MHz. The linewidths necessary 
for these simulations can be found in [1, 2, 36, 43]. 


V. CONCLUSIONS 

We have demonstrated that the 688-nm transition can be 
used to reduce cycle time and increase atom number in ul¬ 
tracold strontium experiments. For the bosonic isotopes, ap¬ 
plying both a 688-nm and a 679-nm laser to the bMOT in¬ 
creases atom number in the metastable reservoir by up to 65 % 
regardless of loading time. If an experiment requires large 
atom number relative to experimental limits, the trap loading 
time can be reduced by a factor of three or better. The en¬ 
hancement is less for ^^Sr due to complications arising from 
hyperfine structure and smaller Lande g-factors for the ^£2 
state. If a second frequency component to simultaneously 
pump |F= 11/2) ^ |F' = 11/2) and |F = 9/2) ^ |F' = ll/2) 
were added to the depumping beam, we believe performance 
comparable to the bosonic isotopes would be achievable. This 
improvement might make the depumping technique a useful 
method to reduce dead time in *^Sr atomic clocks [44]. 

Comparison with a one-dimensional rate equation model 
shows that our results for the bosons are consist with expec¬ 
tations. The initial prediction of 3x increased atom number, 
based on the branching ratio from ^D 2 into ^£2 and ^Pi, is not 
feasible. Simulations of alternative enhancement schemes in¬ 
dicate that pumping on either the 5s5p ^P\ —> 5s5d^D2 tran¬ 
sition or the 5s5p ^Pi —> 5s6d^D2 transition, which are also 
accessible with diode lasers, offers similar performance to the 
approach pursued in this work. Regardless of the exact imple¬ 
mentation, the trap loading enhancement scheme can substan¬ 
tially increase atom number independent of the bMOT loading 
rate or vacuum lifetime. We expect that this method will be 
helpful for experiments benefitting from high atom number or 
faster cycle times. 
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